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Part  I:  Fractional  Dimension  (Fractals,  Bioinspired,  Intelligent  C.) 
sine  wave  O  versus  Weierstrass  Function  C° 


Part  II:  Quantitative  Feedback  Theory 


Part  III:  A  Common  Problem  -  Diffusion  Equation 

(a)  Solve  the  classical  way. 

(b)  Solve  using  Laplace  Transforms. 

(c)  Solve  using  Fractional  Calculus. 

(d)  Examine  Robustness  via  Quantitative  Feedback  Theory. 
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Quantitative  Biofractual  Feedback  Part  I 


.  We  are  now  living  in  a  world  that  is  complex,  distributed, 
but  may  be  highly  vulnerable. 


.  A  better  understanding  of  performance,  and  vulnerability 
of  complex,  distributed  systems  is  required.  How  should 
we  allocate  resources  for  protection? 


The  Part  I  talk  will  have  four  main  components: 

(A)  Pose  the  problem  of  performance  and  vulnerability  in 
complex  and  distributed  networks. 


(B)  Provide  background  material  on  some  pertinent  areas. 


(C)  Using  Computational  Intelligent  methods,  solve  a 
related  problem.  This  will  be  a  “brute  force”  approach. 


(D)  Finally,  hypothesize  some  theoretical  approaches. 


Part  1-A-  Pose  the  Problem: 


Figure  3  -  The  Original  Network-Centric  Distributed  System 


Performance.  Rate  of  flow  through  the  network. 


Vulnerability:  Sensitivity  of  performance  to  attack  of  node. 


(absolute  or 
relative 
objective 
measures 


of  V,  ) 


Part  1-A-  Pose  the  Problem: 

Some  examples  of  important  networks: 

(1)  Power  grids,  railroad  tracks, 

financial  systems,  etc.. 

(2)  Flow  of  people,  water,  food,  medicine. 

(3)  Communication  systems. 

(4)  Information  networks  (Internet), 
email  systems. 

(5)  Physiological  systems  (blood,  oxygen, 
heart  attack,  cell  networks  in  biology). 


(Some  networks  we  may  wish  to  destroy.) 


Part  1-A-  Pose  the  Problem: 


One  Network  we  wish  to  destroy: 


C-IED  Capabilities 


Timeline  for  single  IED  event 


dstl 


21  March  2007 
©  Dstl  2001 


Dstl  is  part  of  the 
Ministry  of  Defence 


A  second  important  network  to  introduce 
congestion  or  denial  of  service: 


Part  1-B-  Background  Material 


Optimization 


Graph 

Theory 


Information 
\  Theory 


Bioinspired 

(Fractals) 
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Part  1-B-  Background  Material 


Bioinsoired  -  Fractals 


Oxygen  Diffusion 
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Part  1-B-  Background  Material 


Bioinspired  -  Fractals 


.  The  Latin  fractus  =  “broken”  or  “fractured” 

.  Fractals  -  scale  free  (self-similar),  irregular  overall  length 
scales,  (self  similar  means  the  structure  is  invariant  to  change  in 
scale).  Forever  continuous  but  nowhere  differentiable. 


.  Fractals  may  have  infinite  circumference  but  finite  area. 


.  Fractals  can  have  finite  volume  and  infinite  area. 


.  A  fractal  can  be  defined  in  the  sense  of  a  recursive  equation: 

Zn+1  =  f(zn) 


.  This  is,  apparently,  the  optimal  way  to  distribute  flow. 


.  Non  Euclidean  Geometry. 

.  Fractal  examples  (trees  (branches),  rivers,  lighting  bolts,  cells, 
lung  passageways,  blood  vessels,  leaf  patterns,  cloud  surfaces, 
molecular  trajectories,  neuron  firing  patterns,  etc.). 


Fractals  -  Lets  Review  the  Area 


Area  1 


B.  Mandelbrot  (1960,s)  asked  the  question:  “How  long  is  the 
coastline  of  Britain?” 


(Suppose  we  measured  the 
coastline  with  a  ruler  that  got 
smaller  and  smaller?) 


A  fractal  has  statistical  self-similarity(  power  law,  self  affine). 

A  fractal  has  N  identical  parts  with  scale  factor  L, 


The  Hausdorff  dimension  is 


Area  =  L2 


Length  =  L 


Volume  =  L3 


(Measurement)  =LD 
implies  log  (Measurement^  D  (log(L)) 


log  (Measurement) 

log  L 


*  Integer 


Fractals  -  Lets  Review  the  Area 


j— IB 


\og(Measurement ) 


(Measurement)  =  L  D 

L  a  A  1/2  a  V1/3 


For  irregular  surfaces,  we  can  define: 

Let  N  =  the  number  of  divisions  of  fixed  length. 
Let  r  =  length  of  a  ruler. 

as  r  — »  0 


D  - 


lo  g{Total  Length) 


\og(\!r) 


Fractals  -  Lets  Review  the  Area 


log  {Measurement) 

log  L 


Total  Length  =  LD  where  1  <  D  <  2 

Koch  Snowflake 


Length  =  4  =  measurement 


Projection  =  topological  dimension  =  3 


igg(4 ; 

log(3) 


1.26185.. 


Fractals  -  Lets  Review  the  Area 
Different  versions  of  the  Koch  snowflake. 


Circumference  |jm  ^0ta|  length)  oo 

=  total  length  n— ^ 


=  (4/3  )n 

Log(total 

length) 


Biofractals 

21  orders  of 
magnitude 

Microbe  =  1CH3  g 

Whale  =  108g 


Power  law 


Log(1/e) 


Fractals  -  Lets  Review  the  Area. 


D  = 


\og(Measurement ) 

logL 


Area  1 


How  to  determine  Measurement? 

We  “cover”  with  boxes  or  disks. 


Fractals  -  Cantor  Set  (Cantor  Dust 


D  = 


\og(Measurement ) 

logL 


Area  1 


(remove  the  middle  third) 


Basic 


2/3 


2/3  (2/3) 


D  =  ]og(2)  =  0  63092 
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What  is  the  Complement  of  the  Cantor  Dust  Set? 

The  set  of  deleted  points  of  Lebesque  measure  1 
The  remaining  points  of  Lebesque  measure  0. 


0.25 


0.50 

Time 


i.oc 


Fractional  Calculus  -  Main  Points 


(non  Euclidean  geometry) 


What  can  n  be? 


(Notation 
invented  by 
Leibniz) 


Answer  (In  1695,  L’Hopital  asked 

Leibniz,  suppose  n=  !4?) 

n  =  integer  =  1,  2,  3  4, 
n  =  negative  integer  =  -1 ,  -2,  -3 
n  can  be  a  non  integer,  n  =  14,  5/6. 
n  can  be  a  negative  non  integer,  n  =  -.6,  -3.4, 


n  can  be  irrational:  fj  — 


n  can  be  a  complex  number: 


n  =  V— 1 


Fractional  Calculus  -  Main  Points 


(non  Euclidean  geometry) 


Why  Study  Fractional  Calculus? 


Composite  Materials 


Log  of  frequency 


CO 


Fractional  Calculus  -Main  Points 
Why  use  Fractional  Calculus? 

(1 )  It  can  deal  with  functions  that  are  forever  continuous 
and  nowhere  differentiable  (fractals). 


(2)  It  has  the  property  of  self  similarity  (scale  invariance) 

d5l2(ccy)  dV2(ay)  dm(ay)  dV2(au)  dU2(au) 
d(at)5'2  d(at)2'2  d(at)V2  d(aifn  d(at)V2 

dqf(bx)  _  d“f{bx) 

[dx]q  [d(bx)]‘> 

(3)  It  is  also  of  the  form: 


zn+1  =  f(zn) 

(Iterated  function  theory). 

(4)  It  can  also  solve  partial  differential  equations: 


d2u(x,t) 

dx2 


du(x,t ) 


dt 


Fractional  Calculus 


An  Easier  Way  to  View  the  Self  Similarity  Property 


A  power  law  f(x)  =  xa  has  the 
property  that  that  the  relative 
change  in  / (fa;) 


/(*) 


=  k 


a 


Is  independent  of  x 


In  this  sense,  the  functions  lacks  characteristic  scale 
(scale  free  or  scale  invariant).  Let  us  evaluate  fO™) 

f  O) 

Let  x  =  ya 
Then 

f(kx) 
fix) 


Note:  no  dependence  on  x 


Fractional  Calculus  -Main  Points 


(310  year  old  area).  Non  Euclidean 

Common  Properties 


(1 )  Scale  Invariance  -  Self  Similarity. 

d“f{bx)  _  d“f{bx) 

[dx]q  [d(bx)Y 

(2)  Weierstrass  Function: 

oo 

/(*) = Z  a"  cos (bn7a) 

n=0 


(3)  Solves  Systems  in  Nature  (Diffusion  equation). 


Fractional  Calculus  -Other  Points 


(310  year  old  area).  Non  Euclidean 


Forever  continuous  nowhere  differentiable. 


Weierstrass  Function: 


/(*) = z  an  cos (bn7uc) 

77 =0 


0<  a  <1,  b  is  a  positive  integer  and  ab  >  1  +  (312)% 


Solves  Systems  in  Nature  (Diffusion  equation). 


Fractional  Calculus  -Other  Points 


Weierstrass  Function  (Why?):  mm 

00 

/(*) = X  an  cos {bn7uc) 

n=0 

0<  a  <1,  b  is  a  positive  integer  and  ab  >  1  +  (3/2)rc 


Step  1 :  We  understand  the  radius  of  convergence: 

1  ,  2  3  4  5  6 

- —  1  +  x  +  x  +  x  +  x  +  x  +  x  +  ... 

1  -  X 

00  J 

^  x  = - =  1  +  x  +  x  +  x  + ... 

n=o  1  -  * 


Fractional  Calculus  -Main  Points 

(Solution  of  the  Diffusion  Equation) 


Area  2 


00 


T(z)  =  [e  Uuz  xdu ,  r(l)  —  1,  r (z  +  1)  —  zT(z), 

0  Thus:  (  r (z  +  1)  =  z!,  )  (  T(— )  = 


Step  1 

d 


-  Derivatives  in  X 


m 


m 


dx 

l 

dxp 


x  -mx 


m — 1 


dp 


xm  = 


m! 


.m 


X  = 


dx 13  ( m  — P)\ 

j. 

r(m  +  l)  oP 


x 


m-f3 


but  (3  may  not  be  an  integer 


T(m  -  p  + 1) 


x 


<ix^ 


X  = 


r(i  + 1)  A  2 


r(i--  +  i) 
'  2 


X  2  = 


X 


This  now  generalizes  for  derivatives  in  eax 


ax 


Dv  eax  =  av  e 

(v  not  an  integer) 

Generalizations  to  functions  that  can  be  written  in  a  power  series: 


q 

/,w=y«, +b„x 


n 


n  =  0 

Generalizations  to  functions  that  can  be  written  in  an  exponential 
series:  „  JO 


q 


f2(t)  =  Yjan  +  V 


n 


e 1  =  cos  {6)  +  i  sin(#) 


n- 0 


id  -id 

e  +e 


Euler’s  Law: 


cos(0) 
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Fractional  Calculus  -Main  Points 

(Solution  of  the  Diffusion  Equation) 


^Areall 


Step  2  -  Laplace  Transform 


GO 


F(s)  =  L[f(t )]  =  J  f(t)e-s'dt 


Then: 

-l 


which  holds  if 


e“  |  /(f)  |<  M  <  oo 


L-[F  «]  =  /«) 


-1 

1  t>  L-\\]  =  — 


1 


r(/?  + 1) 


i  i 

r  r(—)  (a/T)?2 


Step  3  -  Diffusion  Equation: 


d2u(x,t)  2  du(x,t) 


00 


dr 


=  a 


dt 


U (x,  s)  =  L[u(x ,  0]  =  J  ^  **w(x,  t)dt 


TVdu  1  S2^  \  ,  1  52C/  A 

£[— - ^  (x,s)-  / (x) — =  0 


dt  a  dx 


a  dx 


1  1  *  GO 

U(x,s)  =  Aexasi  +  ^  f  e“^x'r|/(r)^r 

a  2js  i_ 


1 


oo  -(x-ry 


u(x,  t )  =  — ^  J  e  4'  /(r) d r 


Part  1-B-  Background  Material 


Information  Theory 


H(x)  = 
Input 

Information 


H(x  /  y)  =  Lost 
Information  = 
Equivocation 


l(x;y)=  Mutual 
Information 


H(y)  = 

Output 

Information 

— ► 


H(y/x)  =  Noise  or 


Spurious  Information 


Dr  -  H(x/y)  +  H(y/x)  (metric  not  a  measure) 

p(x,y)>  0  for  all  x  and  y.  (non  negativity) 
p(x,y)  =  p(y,x)  (symmetry) 

p(x,z)  £_  p(x>y)  +  p(y,z)  (triangular  inequality) 
p(x,y)  =  0  IFF  x=y  (identity  of  indiscernibles) 


Mutual  Information  (l(x,y))  is  well  embraced  by  numerous 
disciplines.  (Ml  is  the  reduction  in  uncertainty  in  an  input 
object  by  observing  an  output  object). 


Part  1-B-  Background  Material 


Information  Theory 


H(x)  = 
Input 


Information 


H(x  /  y)  =  Lost 
Information  = 
Equivocation 


l(x;y)=  Mutual 
Information 


H(y)  = 
Output 

Information 

. ■■■■■■■■»• 


H(y/x)  =  Noise  or 
Spurious  Information 


Why  are  we  interested  in  flow  rate? 

Units  of  l(x;y)  are  bits/sec 
Therefore  bits  =  l(x;y)  At  where 
At  =  time  to  complete  a  task. 

Suppose  we  view  bits  as  discrete  events. 

If  bits  =  events  =  fixed  then: 

min  l(x;y)  =>  max  At,  max  l(x;y)  ^  min  At 

Optimal  Performance  Optimal  Network  A  ttack 


At  = 


events 
I(x;  y ) 


Part  1-B-  Background  Material-Graph  Theor 


(1)  Random  Graphs.  (Less  vulnerable,  uniformly 
connected). 


(2)  Scale  free  graphs.  (Highly  vulnerable,  not 
uniformly  connected). 


»ii 


Part  1-B-  Background  MateriaJ 


Graph  Theory  (Spatial  Construct) 


The  Internet 


Internet-Map 


Part  1-B  -  Background  Material 

Graph  Theory  (Spatial  and  Temporal) 

The  Internet  is  dynamically  scale  free  (evidence) : 

IEEE/ACM  TRANSACTIONS  ON  NETWORKING,  VOL.  2,  NO.  I,  FEBRUARY  1994 


AT  =  100  sec. 

^  B0OOO 
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Reference:  W.  E.  Leland,  et  al.,  lEEE/ACM  Trans,  on  Networking,  vol  2,  no.  1, 
Feb.  1994,  “On  the  Self-Similar  Nature  of  Ethernet  Traffic  (Extended  Version).” 


Part  1-B  -  Background  Material 


Graph  Theory 


The  Internet  is  dynamically  scale  free  (evidence) : 


JRQ VELLA  AND  BESTAVRQS:  SELF-SIMILARITY  IN  WWW  TRAFFIC 
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R0f0j-0PQ0;  M.  Crovella  and  A.  Bestavous,  “Self-Similarity  in  World  Wide  Web 
Traffic:  Evidence  and  Possible  Causes,”  IEEE/ACM  Trans.  On 
Networking,  Vol.  5,  no.  6,  December,  1997. 


Other  Physiological  Evidence 


Heartbeat  Intervals 


Figure  2.  The  time  series  of  heartbeat  intervals  of  a  healthy  young  adult  male  is  shown.  It  is  clear 
that  the  variation  in  the  time  interval  between  beats  is  relatively  modest,  but  certainly  not  negligible. 


Heartbeat 

intervals 


B.  J.  West,  “Fractal  Physiology,  Complexity,  and  the  Fractional  Calculus,” 
Chapter  6,  in  “Fractals,  Diffusion  and  Relaxation  in  Disordered  Complex 
Reference!  Systems,”  in  Advances  in  Chemical  Physics,  vol  133,  part  B,  John  Wiley,  2006, 
Eds.  W.  T.  Coffey  and  Y.  P.  Kalmykov. 


Other  Physiological  Evidence 


\CIJ)  M.VJ,  ailU  L  — 7  w,  U1C  ICUglll  Ul  UlC  J/nawM.  mv  aM66v««. 

Cantor  set  exponentially  goes  to  zero,  the  fluctuations  of  the  data.  This  prop 


A.  NORMAL  SINUS  RHYTHM  TIME  SERIES 


73  1-0  I 


0.5 


4.  Computer  simulation  of  a  fractal  lung,  in  which  the  boundary  conditions 
influence  morphogenesis.  The  boundary  was  derived  from  a  chest  radiograph.  The 
model  data  are  in  good  agreement  with  actual  structural  data  191. 


I 


Sinus  Rhythm 
Intervals 


B  NORMAL  SINUS  RHYRHM: 

INVERSE  POWER-LAW  SPEC TRUM 


Reference: 

W.  Deering  and  B.  J.  West,  “Fractal  Physiology,”  IEEE  Eng.  In 
Medicine  and  Biology,  June,  1992. 


Additional  Background  Material  -  H.  Jeong  -  Complex  ‘07 
(The  difference  between  random  and  scale-free  graphs) 


Highway  network 
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Airline  network 


"Number  of  nodes  with  k  links 


Mathematically?  via  Degree  distribution  P(k) 


Random  Scale  Free 


Bell  Curve 
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Power  Law  Distribution 


World  Wide  Web 

Node(point):  web-page  ^ 
link(line):  hyper-link 


The  problem  is  to  discern  (for  each  application): 


(1  )What  are  the  nodes? 
(2)What  are  the  links? 


INTERNET  BACKBONE 

Nodes:  computers,  routers 

Links:  physical  lines 


(Faloutsos,  Faloutsos  and  Faloutsos,  1999) 


Cumulative  distribution,  P(k 


SEX-Web 

Nodes:  people  (females;  males) 
Links:  sexual  relationships 


(Liljeros  et  al.  Nature  2001) 


wtjsum® 
mn  =m 


O  Male  (Area4) 

O  Female 


Sexual  Relationships  in  Jefferson  High  School 


ACTOR  CONNECTIVITIES 


Nodes:  actors 


Links:  cast  jointly 


IMP-  Internet  Movie  Database 


requiehO 


Days  of  Thunder  ( 1990) 
Far  and  Away  ( 1992) 
Eyes  Wide  Shut  (1999) 


N  =  212,250  actors  <k>  =  28.78 


P(k)  ~k-T  y=2.3 


IO'1 


io'2 

10'3 

io'4 

10'5 

io'6 


10°  io1 


io2  io3 


k 


SCIENCE  CITATION  INDEX 


Nodes:  papers 
Links:  citations 

1736  PRL  papers  (1988) 


P(k)  ~k^ 

(y  =  3) 

(s. 


v\ 

\  \ 

*  *  *  ?4  7 
1  /// 

KPZE 
PRL  56,  8 

/  /  / 

quation 

189  (1986) 

m, 

m 

m 

M 

12  34 . 1247 

Redner,  1998) 
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SCIENCE  COAUTHORSHIP 


(collaboration  network) 


Nodes:  scientist  (authors) 
Links:  write  paper  together 
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Other  Examples  of  Scale-Free  Networks 

Email  network 

Nodes:  individual  email  address 
Links:  email  communication 

Phone-call  networks 

Nodes:  phone-number 
Links:  completed  phone  call 

(Abello  et  al,  1999) 


Networks  in  linguistics 

Nodes:  words 

Links:  appear  next  or  one  word  apart  from  each  other 

(Ferrer  et  al,  2001) 

Networks  in  Electronic  auction  (eBay) 

Nodes:  agents,  individuals 
Links:  bids  for  the  same  item 


(H.  Jeong  et  al,  2001) 


THEN  WHY?? 


Area  4 


(i)  Efficiency  of  resource  usage. 

Diameter  (Scale-free)  <  Diameter  (Exponential) 

(*  Diameter  ~  average  path  length  between  two  nodes) 


(ii)  Robustness  of  complex  networks. 

Scale-free  networks  are  more  robust  under  random 
errors,  but  very  vulnerable  under  intentional  attacks! 


Scale-free  Networks  are  efficient/robust. 


Points: 

(1  Vulnerability  ([robustness]'1)  is 
predicated  on: 

(a)  Architecture  of  network 

(b)  Type  of  attack. 


What  is  the  Real  Problem? 


Area  4 


Most  networks  are  not  static,  they’re  dynamic! 


e.g.  real  metabolic 
networks  are  DYNAMIC!! 


Let  us  stop  with 

Graph  Theory  and 
move  on  to  the  last 
area 


Optimization  Theory. 


Part  1  -C  -  Let  us  work  a  practical 
example: 

(15  flows) 


Structure  For  the  CAPS  Simulation  using  GAs 


Minimum  (8  links) 


Maximum  (20  links) 


Case  1  -  Only 
Communications 
Through  the  ATOF 


Case  2  -  Full  Spam 
Communications 


Part  1-C  -  Issues  of  Vulnerability  and  Performance 


KirchhofFs  Law  and  Cut  sets 


I  Currents  =  0  into  a  node. 


Kirchoff’s  Law  also  applies  in  Graph  Theory 


Part  1  -C  —  Issues  of  Vulnerability  and  Performance 


Kirchhoff's  Law  and  Cut  sets 


Maximum  Flow 


Minimal  Flow 


Cut  set:  flows  in  =  flows  out 
=  10  units 


Cut  set:  flows  in  =  flows  out 
=  1  unit 


ATOF: 

fx+  f2+  f4+  f6+f8  =  f1  +  f3+  h  +  f7 

PS: 

fn+  f1  =f9  +  f2  +  f12 

RS: 

f15+  f7+f9+f13=  f10+f8 

FS: 

'10  +  '12+  * 3  =  ' 14+  *13+  <4+  'll 

CS: 

'5  +  '14  =  '15+  '6 

AT 

Sensitivity  = 

t  r  dT  W 

Sw  :=  lim  J  =  (7"  ^0) 

AW—>o  AW  dW  T 

w 

Let  T  =  cut  set  flow,  let  W be  the  Ml  =  l(x;y). 


Area  5 


j  =  1,  11  free  chromosomes 

3  bit  word  for  each  chromosome. 


j  =  i 

j  =  2 


0 

1 

1 

1 

1 

0 

0 

1 

j  =  11 


(811  possibilities,  NP  Hard) 

Fig.  9  Configuration  for  the  Chromosome 


How  the  Optimization  is  Conducted  (Elite  Pool) 


Bump  Out 


Bump  In 


Fig.  10  -  Maximizing  (l(x;y))  vs.  Pool  Entrance  Number 

1.4 
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1.1 
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0.9 
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0  50  100  150  200  250 

Fig.  11  -  l(x;y)  Minimization  vs.  Pool  Entrance  Number 


Sensitivity  Results  -  Logistics  Problem 


Area  5 


Sensitivity  function  in  equation  (32)  for  5  computer  runs 


ATOF  vs  PS  for  5  computer  simulation  runs 


Figure  (12)  -The  sensitivity  Function  defined  in  equation  (32)  for  ATOF  vs  PS 


Simulation  is 
Sometimes  termed 


“Experimental 
Mathematics ” 


Other  Common  Intersections 

Causality  Map 


Part  D  -What  is  the  solution  in  a  theoretical  sense? 


.  Bioinspired  =>  Perhaps  we  should  not  think  Euclidean? 

.  Fractional  Calculus  may  capture  dynamics. 

.  Here  may  be  a  hypothesized  solution? 

Network  Science 


Robotics 


Minimize  (J.,) 


Figure  3  -  The  Original  Network-Centric  Distributed  System 

Minimize/Maximize  (l(x;y)) 


Subject  to  constraints: 

x  =  je 


Subject  to  Constraints: 

m=  o 


J5/2(oy)  J3/2(oy)  J1/2(cy)  _  d^jau)  d^jau) 

d(at)512  d(atyi2  d(at)V2  d(at)312  d(at)112 


End  of  Part  I  -  Quantitative  Biofractal  Feedback 


.  Performance  and  vulnerability  of  distributed 
systems  needs  to  be  objectively  quantified. 


.  We  can  learn  from  biological  systems  (fractals). 
Also  the  fractional  calculus  may  offer  a  venue 
to  characterize  dynamics. 


.  There  are  many  common  connections 
between  five  different  areas.  For  example, 
the  diffusion  equation  is  bioinspired. 


.  Computational  methods  allow  us  to 
synthesize  a  brute  force  approach  for  insight. 


.  Much  more  work  needs  to  be  accomplished. 


Part  II  -  Brief  Review  of  QFT 

.  Quantitative  Feedback  Theory  originated  in  the  1960’s  by  Isaac 
Horowitz  using  frequency  domain  methods  for  efficient  robust 
control  design.  In  1972  a  seminal  paper  was  published. 

.  QFT  has  been  used  in  Flight  Control,  Robotics,  Power  Systems, 
unmanned  air  vehicles,  and  many  other  applications. 

.  The  controller  is  determined  by  a  loop  shaping  process  employing 
a  Nichols’  Chart  that  displays  the  stability,  performance  and 
disturbance  rejection  bands. 

.  A  typical  QFT  Controller  (synthesis)  satisfies  certain  attributes: 

(a)  Robust  Stability. 

(b)  Reference  Tracking. 

(c)  Disturbance  Rejection. 


1  DoF  System 

D, 


QFT  Basics 


In  the  Absence  of  Disturbances  D|  and  D0: 

Let:  L  =  Loop  Gain:  L  —  C  P 

Then  the  closed  loop  transfer  function  between  Y  and  R  is: 

TQ0  =  T,,  =  F(s)L(s)  =  Output^ 
R(s)  1  +  L(s )  Input 


The  Sensitivity  of  The  Closed  Loop  Transfer  Function  T(s) 
to  plant  variations  P(s)  can  be  specified  via: 


1 

1  +  L(s) 


RM 


QFT  Basics 

n.  n 


For  QFT  Design,  we  have  at  least  3  criteria  to  meet: 

(1)  Robust  Stability  (closed  loop  Robust  Stability) 


m 

1  +  L(s) 


=>  This  is  a  constraint  on  the  peak  magnitude  of  the 
closed  loop  frequency  response. 

(2)  Reference  Tracking.  Let  TLand  Ty  be  the  upper  and 
lower  transfer  functions,  then  we  require: 


|TL(jco)|  <  |  TO  co)  |  <  T|j(jco)| 


(3)  Disturbance  Rejection:  We  require: 


1 

1  +  L(jco) 


1 

"  W(jco) 


Where  W(jco)  is  a  weighting  function  (of  frequency). 


Note  conditions  (1-3)  are  for  the  class  of  plants  P  s  {P,} 


RM 


QFT  Basics 

n.  n 


For  the  Disturbances  Dj  and  D0 


The  Transfer  Function  between  Djand  Y  is  given  by: 


Y(j<s)  PUto) 
D,  ( jco )  1  +  L(jco) 


The  Transfer  Function  between  D0  and  Y  is  given  by: 

T  YU  to)  1 

D0  a  to)  1  +  L(jto) 

Then  the  Disturbance  Rejection  Can  Be  Specified  via: 


T,  <B 


di 


di 


Where  the  Bdj  and  Bd0  are  frequency  dependent  functions. 
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Applications 


Part  III  -  The  Diffusion  Equation 

Why? 

(1)  Many  biological  systems  can  be  characterized  in  this  manner. 

(2)  Outside  biology,  diffusion  is  a  fundamental  process  (thermal 
chemical,  other  physical  processes  of  all  types). 

(3)  The  diffusion  equation  satisfies  a  fractional  differential 
equation. 

(4)  The  diffusion  equation  is  also  a  type  of  fractal. 

Consider  the  following  physical  problem: 

Let  u(x,t)  be  the  temperature  distribution  in  a  cylindrical  bar  of 
finite  length  L  oriented  along  the  x-axis  and  perfectly  insulated 


and  k  is  the  thermal  conductivity,  c  is  the  specific  heat  and  5  is 
the  linear  density  (mass/unit  length). 


The  initial  condition  is:  u(x,0)  =  f(x) 

The  boundary  conditions  are:  u(0,t)  =  0  =  u(L,t)  V  t 


Part  III  -  The  Diffusion  Equation 

d2u(x,t ) 

dx2  dt 

u(0,t)  =  0  =  u(L,t)  Boundary  Conditions 
u(x,0)  =  f(x)  Initial  Condition 

Possible  ways  to  solve  the  equation: 

(1)  Fourier  Method  -  Separation  of  Variables. 

(2)  Laplace  Transforms. 

(3)  Fractional  Calculus. 


Now  examine  Robustness  via  Quantitative  Feedback  Theory 


Part  III  -  The  Diffusion  Equation 


d2u(x,t )  _  2  du(x,t ) 

A  CL 


Initial  Condition  u(x,0)  -  f(x) 


dx 


dt 


Boundary  Conditions:  u(0,t)  =  0  =  u(L,t) 

(1)  Fourier  Method  -  Separation  of  Variables. 

Assu  me  u(x,t)  =  X(x)T(t) 


a2X(x)T(t)  =  T(t)X"(x ) 


f!M  =  Z^)=  constant  =  -2 

Tit)  X(x) 


T(t)  =  -(A/a2)T(t ) 


T{t)  =  Ae 


—  )( 


and  X"(x)  =  -AXix)  ^ 

=>  X(x)  =  Bsm(JIx)  +  C^NXx)  but  u(0,t  )=  0=>  C=0 


uAx,t)  =  T  it)XAx) 


and  u{x,t)  =  2>  ,  (x,  t) 


00 


(X,  t)  =  X  A,(sin(^)  (a2L2) ) 


( n27i2t ) 


u 


11— \ 


L 


Dn  =  —  [  /(.x)  sin 
L 


^  htdC 


0 


V 


L 


dx 


J 


(Can  show  the  infinite  series  converges) 


Part  III  -  The  Diffusion  Equation 


d2u(x,t )  du(x,t ) 


dx 


dt 


Initial  Condition  u(x,0)  -  f(x) 


u(x,t)  bounded,  t>  0,  -oo  <  x  <  oo 
(2)  Laplace  Transforms. 


OO 


Define  the  Laplace  Transform  Variable:  U(x,s )  =  \e  tsu{x,t)dt 

oo  ^  J 

-ts  ou 


J -ftdt  =  sU(X'S}  ~  M(x’0)  =  s  U(x,s)  - 


f(x) 


0  , 

du  d  u 

If  —  and  — r  are  bounded  and  continuous 


dx 


dxz 


uu 

J 


„  dlu  ,  dlU 
e  st  — xdt  = 


dxz 


dx ' 


Now  Laplace  transform  the  partial  differential  equation 


o  =  1 


du  82u 


dt  dx2 

and  solve  for  U(x,s)\ 


=  sU(x9s)  —  f(x)  - 


d2U 


1 

U  (x,s)  =  f  e~rAx~y\ f  {y)dy 

2  4s  L 


Forcing 

Function 


To  find  u(x,t),  we  need  to  find  the  inverse  Laplace  transform 

u(x,t)  =  L_1[  U(x,s)  ] 

or  T  i  +0° 

J  e-^-yl  f(y)dy 


u(x,t)=  L-i 


By  integration  in  the  complex  plane  we  can  show: 


u(x,t)  =  , -  i  e 

V2  7Tt 


*  uu 

1  J 


fiyVfy 


—  00 


Part  III  -  The  Diffusion  Equation 


d2u(x,t )  _  2  du(x,t ) 

a  ^ 


dx‘ 


dt 


Boundary  Conditions: 


Initial  Condition  u(x,0)  -  f(x) 


u(0,t)  =  0  =  u(L,t) 


(Heaviside  Operational  Calculus) 

d2u  9  du 


Consider: 


=  a 


dx2  dt 

d  82u  2 

a  *  ^  =  apu 

(treat  p  as  a  constant  and  solve  for  x) 


Let  p  = 


u(x,t)  =  Ae 


■apl/2x 


+  Be 


apl/2x 


On  physical  grounds,  B  =  0 

_  -axp 


=>  u]  (x,  f)  =  e 

^  00  / 

Or:  u{x,t)  =  u0  +  ^ 


1/2 


1/ 


0 


(-ax)”  /2 

,  P  uo 


n= 1 


ft! 


(can  ignore  positive  integral  powers  of  p) 


u(x,t)  =  w0 


Part  III  -  The  Diffusion  Equation 


d2u(x,t )  _  2  du(x,t ) 

a  ^ 


dx‘ 


dt 


Initial  Condition  u(x,0)  -  f(x) 


Boundary  Conditions:  u(0,t)  =  0  =  u(L,t) 

Now  examine  Robustness  via  Quantitative  Feedback  Theory 


Step  1 :  Let  us  examine  a  heat  control  problem. 
(Define  units  of  all  quantities  to  generalize. ) 


Step  2:  Let  us  build  a  controller  within  a  QFT  context. 


Step  3:  We  have  now  solved  a  heat  control  problem.  Now 
generalize  to  flow  problems  as  in  networks. 

Again  look  at  the  units  of  all  variables. 


Part  III 

Step  1 :  Let  us  examine  a  heat  control  problem. 

Let  udesired(x,t)  =  desired  temperature  =  uD(t)  (assume  x=const). 
Let  uactua)(x,t)  =  actual  temperature  =  ua(t) 

Temperature  error  eT(t)  =  uD(t)-ua(t) 


Units  Analysis:  U;(t)  =  temperature  -  C° 

C  =  Thermal  Capacitance  =  kilo  cal  /  C° 
q(t)  =  heat  input  -  kilo  cal  /  second 
Rx  =  Thermal  Resistance  -  C°sec  /  kilo  cal 


du 


C  n 

,  —  Qi  ~  Vo  Where:  — 


U 


a 


dt 


^  du 

c — -  +  q0  =  qi 


R 


T 


dt 


du„  u„ 

+—  =  q, 


C""a 


rtc 


dt  Rj 

du 


dt 


+  ua  -  Rrq, 


UM  rt 

6iM  i  +  ktCs 


Part  III 


Step  2:  Let  us  build  a  controller  within  a  QFT  context 


QFT  Goals: 

(1)  Stability  7^) 


R 


■±Q 


L 


Y 


1 +  L 


is  stable. 


■*  Y 


L  =  GP 


(2)  Tracking  Specifications 

l~rLGoo)|  £  |F(j  (D)  T(j  co)|  <  |Tu(j  oo)|  =^>  use  F  for  prefilter. 


(3)  Disturbance  Rejection 
max  |TD(jco)|  <  |MD(jco)| 

QFT  Design  Procedure: 


T  = 

1  Di 


T  = 

1D  0 


p 


1 +  L 
1 

1  +  1 


(a)  Find  the  plant  templates  Ps  {PJ  -  Nichols  chart. 

(b)  Generate  Performance  Bounds  from  Nichols  chart. 

L0(s)  =  P0(s)  G(s) 

(  c  )  Loop  Shaping:  Add  poles  and  zeros  to  L0(s). 

(d)  Design  Prefilter  F  (  keep  |TL|<|  F  T  |  <  ) 

(e)  Finally  to  determine  the  final  controller 


4M 

Pis) 


Done! 


Part  III 


Step  3:  We  have  now  solved  a  heat  control  problem.  Now 
Generalize  to  flow  problems  as  in  networks. 


Heat  Control  Problem: 


The  Network  Flow  Problem 


System  jn?  + 

- ► 


Controller  out 


> 


controller 


? 


Plant  out 


? 


Let  us  review  the  units  of  variables  of  interest: 


Heat  Control  Problem 

u  |  units  of  (C° ) 

q  units  of  (kilo  cal/sec) 
C  units  of  (kilo  cal  /  C°) 


c 


du 


dt 


a  =  d, 


do 


Network  Flow 

System  in  ? 


Controller  out  ? 


Plant  out  ? 


Part  III 


Step  3:  We  have  now  solved  a  heat  control  problem.  Now 
Generalize  to  flow  problems  as  in  networks. 


Suggestions: 

Heat  Control  Problem  -  flow 

q  *  time  =  kilo  calories 


Network  Problem  -  flow 

bits/  sec  *  seconds  =  bits 

events/second  *  seconds 
events 


Equate  the  above  variables  (Ml=q,  events  =  kilo  calories) 


events  - 


bits  =  j  (mutual  information) 


Heat  Control  Problem 

u  |  units  of  (C° ) 

q  units  of  (kilo  cal/sec) 

C  units  of  (kilo  cal  /  C°) 

^  dll 

C-±  =  qt-q0 

dt 


Network  Flow 

System  in  ?  =  J  MI 

Controller  out  7=  ftfj 


Plant 


(Recall  we  modulated  Ml  in  the  example) 
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Part  III 


System  in? 


Controller  out 


> 


controller 


Events/bitsD 


Ml  flow 


Events/bitsA 


> 


controller 


> 


Network  Flow 

System  in  ?  =  f  MI 


Controller  out  ?  —  ][fj 


Plant  out  ? 


(Recall  we  modulated  Ml  in  the  example) 


Summary  and  Conclusions 
Part  I  -  Fractional  Dimensions  - 
non  Euclidean  World. 


Part  II  -  Quantitative  Feedback 
Theory. 


Part  III  -  Diffusion  Equation. 


The  Future  -  Modeling  networks  as 
control  systems  and  applying  these 
techniques.  QFT  helps  because  it 
can  view  robust  control  in  terms  of 
simple  Bode/Nichols  plots. 


